REPORT  DOCUMENTATION  PAGE 

Form  Approved 

OMB  No.  0704-0188 

The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 

gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  the  burden,  to  the  Department  of  Defense,  Executive  Services  and  Communications  Directorate  (0704-0188)  Respondents  should  be  aware 
that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OIS/IB 
control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ORGANIZATION. 

1.  REPORT  DATE  (DD-MM-YYYY)  2.  REPORT  TYPE 

04-02-2008  Conference  Proceeding 

3.  DATES  COVERED  ( From  -  To) 

4.  TITLE  AND  SUBTITLE 

A  Practical  Point  Spread  Model  for  Ocean  Waters 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

0602435N 

6.  AUTHOR(S) 

Weilin  Hou,  Deric  Gray,  Alan  D.  Weidemann,  Robert  A.  Arnone 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

73-6867-07-5 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory 

Oceanography  Division 

Stennis  Space  Center,  MS  39529-5004 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

NRL/PP/73  30-07-72 10 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Office  of  Naval  Research 

800  N.  Quincy  St. 

Arlington,  VA  22217-5660 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

ONR 

11.  SPONSOR/MONITOR’S  REPORT 

NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release,  distribution  is  unlimited. 


13.  SUPPLEMENTARY  NOTES 


14.  ABSTRACT 

The  scattering  properties  of  the  medium  ultimately  determine  the  outcome  of  the  image  transmission.  For  ocean  optics  research,  the 
scattering  properties  are  often  conveniently  described  and  measured  in  general  by  the  scattering  coefficient  (b),  which  determines 
the  possibility  of  a  photon  to  be  scattered  away  from  its  original  traveling  direction  per  unit  length  by  the  medium  molecules, 
constituents  within  (i.e.  particles  [1]),  and  turbulence  [2].  As  we  know,  this  parameter  (b)  is  an  integration  of  the  volume  scattering 
or  phase  function,  B,  which  details  such  probability  by  the  relative  directions  of  incoming  and  out-going  photons  [1].  These  inherent 
optical  properties  (IOP),  although  measured  frequently  due  to  their  important  applications  in  ocean  optics,  especially  in  remote 
sensing,  cannot  be  applied  to  underwater  imaging  issues  directly,  since  they  inherently  reflect  the  chance  of  the  single  scattering. 


15.  SUBJECT  TERMS 

parameters,  constant,  inherent  optical  properties,  PSF 


16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

ABSTRACT 

18.  NUMBER 

nc 

19a.  NAME  OF  RESPONSIBLE  PERSON 

Weilin  Hou 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

Ur 

PAGES 

5 

Unclassified 

Unclassified 

Unclassified 

UL 

19b.  TELEPHONE  NUMBER  (Include  area  code) 

228-688-5257 

Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39  18 


\ 


IV  International  Conference 


Current  Problems  in  Optics 
of  Natural  Waters 


A  PRACTICAL  POINT  SPREAD  MODEL  FOR  OCEAN  WATERS 

W.  Hou,  D.  Gray,  A.  Weidemann,  and  R.  Arnone 

Naval  Research  Laboratory,  Stennis  Space  Center,  MS  39529,  USA 
hou@nrlssc.navy'.  mil 


1.  Background 

The  scattering  properties  of  the  medium  ultimately  determine  the  outcome  of  the  image  transmis¬ 
sion.  For  ocean  optics  research,  the  scattering  properties  are  often  conveniently  described  and  meas¬ 
ured  in  general  by  the  scattering  coefficient  ( b ),  which  determines  the  possibility  of  a  photon  to  be 
scattered  away  from  its  original  traveling  direction  per  unit  length  by  the  medium  molecules,  constitu¬ 
ents  within  (i.e.  particles  [1]),  and  turbulence  [2].  As  we  know,  this  parameter  ( b )  is  an  integration  of 
the  volume  scattering  or  phase  function,  p,  which  details  such  probability  by  the  relative  directions  of 
incoming  and  out-going  photons  [1].  These  inherent  optical  properties  (IOP),  although  measured  fre¬ 
quently  due  to  their  important  applications  in  ocean  optics,  especially  in  remote  sensing,  cannot  be 

applied  to  underwater  imaging  issues  directly,  since  they  inherently  reflect  the  chance  of  the  single 
scattering. 

The  point  spread  function  (PSF)  describes  the  system  response  to  a  point  source  in  the  medium, 
which  includes  the  effect  of  multiple  scattering.  It  is  the  ideal  parameter  to  study  image  transmission, 
optical  sounding  [3],  reversion  of  transmission  effects,  and  retrieval  of  optical  properties  [4],  Gener¬ 
ally  speaking,  a  2-dimentional  image  of  an  object  is  the  combination  of  original  signal,  f(x,y),  con¬ 
volved  by  the  imaging  system  response  of  a  point  source,  the  point  spread  function  or  PSF  h(x  y),  inte¬ 
grated  over  sensor  space  E: 

£(*>  y)  =  ,  y, )h(x  -xi,y-yi )dx,dyi  .  ( l ) 

The  system  response  includes  those  from  both  the  imaging  system  itself,  as  well  as  the  effects  of 
the  medium  (water  in  our  case).  With  known  characteristics  of  the  imaging  system  and  correct  model¬ 
ing  of the  medium,  theoretically  it  is  possible  to  recover  the  original  signal  by  reversion  or  deconvolu¬ 
tion  [4  5].  Mathematically  PSF  is  equivalent  to  the  beam  spread  function  (BSF)  [6]  which  can  be  eas¬ 
ier  modeled  and  measured. 

Duntley  [7]  reported  extensive  lab  measurements  of  BSFs  of  simulated  ocean  waters  with  re¬ 
markably  different  optical  lengths  (0.5  to  21),  and  summarized  their  findings  in  a  single,  albeit  com¬ 
plex  empirical  relationship: 
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£  =  (13.75- 0.501$-)  -  (0.626  -  0.357^)r  +  (0.01 258  +  0.00354$-)r2 , 
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Notice  in  above  equations  the  optical  length  defined  as  r  =  cr  (c  being  the  total  attenuation  coeffi¬ 
cient),  while  the  single  scattering  albedo  is  6>0  =  b/c  .  $•  =  1/(1  -  a>0)  is  used  to  simplify  the  form. 

Another  empirical  effort  was  carried  out  by  Voss  [8],  fitting  field  measured  data  from  3  different 
types  of  oceanic  waters  (TOTO,  Pacific  and  Sargasso  Sea)  to  the  following  form: 


PSF {6)  =  Bx0~m ,  (3) 

where  Bx  is  a  constant  and  is  the  slope  of  log(PSF)  vs  log(5).  While  m  is  not  a  constant,  but  rather 
a  function  of  IOP  and  r,  the  above  formula  can  fit  all  data  to  less  than  15%  error.  The  m  values  ranged 
from  0.4  to  2.0  as  a  function  of  r(0  to  10).  Other  than  graphical  results,  there  is  no  definitive  relation¬ 
ship  provided  to  allow  comparison  of  this  method  to  other  approaches.  Nonetheless,  this  is  encourag- 


ing  result,  since  such  simple  relationship  can  be  rather  beneficial  to  imaging  needs,  especially  when 
high  frame  rate,  per  pixel  calculations  are  desired  and  hardware  implementation  is  needed.  The  aim  of 
paper  is  to  examine  the  relationship  of  the  PSF  to  commonly  measured  IOPs  in  search  of  a  simpler 
solution  for  real-time  imaging  needs,  rather  than  thoroughly  review  PSF  models  (such  as  some  of  the 
more  sophisticated  analytical  models  [9-11]),  although  comparisons  will  be  made  amongst  measured 
PSFs  to  validate  our  results. 


2.  Theory 

The  above  mentioned  empirical  results,  while  useful,  lacks  the  benefit  of  physical  interpretations 
of  each  parameters  involved.  Analytical  solutions,  while  at  times  complex  or  approximated,  offer  clear 
answers  to  physics  involved  and  help  to  check  numerical  solutions. 

Without  including  all  the  details  here,  Wells  [12]  demonstrated  that  under  small  angle  scattering 
approximation  (SAA),  and  a  phase  function  in  the  following  form: 


R(P\  —  K 

PK  J  2n(9\  +  (?2)3/2  ’ 

the  modulation  transfer  function  (MTF)  can  be  expressed  in  a  closed-form  as 

(4) 

=  e-Wr,'> 
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where  6q  is  related  to  the  mean  square  angle  (MSA).  The  PSF  can  be  calculated  numerically  using 
above  relationships  via  inverse  Fourier  transform,  or  analytically  can  be  expressed  in  1st  order  terms  as 
those  formulated  by  Alan  Gordon  [12] 
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From  (4)  and  (9),  we  approximate  the  integration  with  the  following  form: 

PSF (9)  =  K  hre  '  = 

2  n9n  2  n9n 

(10) 

where  K  is  a  constant  and  does  not  affect  imaging  needs  since  relative  units  of  PSF  are  used.  By  com¬ 
paring  to  simple  numerical  calculations  using  Fourier  transformed  MTF,  we  found  that  n  =  1/(Oo-t0o 
gives  the  best  result  over  different  parameter  ranges.  We  will  further  compare  these  with  empirical  and 
measured  results. 


3.  Results 

Data  were  obtained  during  an  April-May  2006  NATO  trial  experiment  in  Panama  City,  Florida. 
The  amount  of  scattering  and  absorption  were  controlled  by  introducing  Maalox  and  absorption  dye 
respectively  into  a  tap-water  filled  pool.  In-water  optical  properties  during  the  experiment  were  meas¬ 
ured.  These  included  the  absorption’ and  attenuation  coefficients  (Wetlabs  ac-9),  particle  size  distribu¬ 
tions  (Sequoia  Scientific  LISST-100),  and  volume  scattering  functions  (multi-spectral  volume  scatter¬ 
ing  meter  or  MVSM). 

The  measured  volume  scattering  functions  are  used  in  Monte  Carlo  (MC)  simulations  to  derive 
the  PSFs,  that  are  used  in  comparison  to  other  approaches.  The  main  computer  code  uses  MC  to  calcu¬ 
late  the  PSF  using  only  the  inherent  optical  properties.  This  method  is  essentially  a  numerical  experi¬ 
ment  of  the  physical  process  under  consideration,  and  is  based  entirely  on  first  principle  calculations 
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which  divorce  it  from  any  other  assumptions  or  theories.  The  disadvantage  of  this  generality  is  the 
long  computer  time  required,  since  the  results  are  statistical  in  nature  and  each  set  of  parameters  must 
be  run  separately.  The  MC  method  has  in  fact  been  used  previously  to  test  the  range  dependence  of  the 
Wells  theory  [13].  Now  rather  than  calculate  the  PSF  directly,  the  mathematically  equivalent  but  com¬ 
putationally  simpler  beam  spread  function  (BSF)  is  calculated  instead.  The  BSF  is  defined  as  the  nor¬ 
malized  irradiance  distribution  on  a  spherical  surface  of  radius  R  centered  on  a  unidirectional  trans¬ 
mitter,  and  this  definition  provides  a  prescription  for  the  Monte  Carlo  method  to  follow.  Given  the 
IOPs  (phase  function,  absorption  and  scattering  coefficients),  the  radius  of  the  sphere  R ,  and  an  angu¬ 
lar  resolution  A6 ,  the  code  produces  the  irradiance  for  each  angular  bin  over  the  entire  sphere. 

Comparison  of  PSFs  from  different  methods  mentioned  (i.e.,  current  study,  MC  simulated  PSF  from 
measured  MVSM  phase  functions,  and  Duntley’s)  are  shown  in  Figs.  1,  2,  and  parameters  listed  in  Table  1. 


fr0™  di,fferent  methods:  “Thi5”  (current  study),  “Duntley”  (empirical  relationship 

maS}‘  x  (Monte  Car,°  sunulated  PSF  based  on  measured  volume  scattering  function)  during 

NATO  expenment.  Normalized  to  1 00  mrad.  See  text  for  details. 
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Fig.  2.  An  example  of  PSF  comparison  under  different  0O  values.  Notice  that  a  higher  value  of  the  mean  square 
angle  value  seems  to  work  better.  This  is  true  for  other  situations  (not  shown)  as  well.  Normalized  to  100  mrad. 


Table  1.  Optical  parameters  used  to  compare  different  PSF  models  as  shown  in  Figs.  1-2 


Date 

optical 
length  (x) 

measured 
c  (nf ') 

scattering  albedo 

((0o) 

4/27/2007 

3.8 

0.78 

0.87 

4/27/2007 

5.5 

0.78 

0.87 

4/30/2007 

0.88 

0.22 

0.67 

5/01/2007 

7.3 

1.38 

0.89 

5/2/2007 

2.5 

0.91 

0.64 

5/2/2007 

5.5 

0.91 

0.64 

The  current  approach  matches  well  with  both  measurements  under  different  testing  conditions 
shown,  and  can  be  considered  an  able  substitute  for  imaging  processing  purposes  [13].  One  may  notice 
from  Figs.  1,  2  differences  can  be  observed  between  the  two  measured  results,  and  that  of  MVSM  is 
always  higher  than  Duntley’s.  This  is  likely  due  to  the  difference  in  scattering  agents  involved,  as  the 
particle  sizes  will  be  different  by  the  different  preparation  process  (recall  that  Duntley  used  high  purity 
filtered  water  as  their  experiment  base  in  the  tank,  before  adding  similar  absorption  and  scattering 
agents).  In  fact,  the  differences  in  PSF  can  be  as  high  as  an  order  of  magnitude  under  certain  condi¬ 
tions  at  small  angles  (eg.  coo=  0.87,  x  =  5.5),  or  ±1  in  PSF  log  slope.  Examination  of  the  parameter  re¬ 
lated  to  0O  suggest  that  a  higher  value  (0.1)  is  more  appropriate  for  the  data  involved  (Fig.  2).  Initial 
testing  with  imagery  obtained  during  NATO  experiment  show  little  differences  using  any  of  the  PSFs 
under  typical  conditions  tested  [4,  14].  It  is  worth  pointing  out,  however,  that  the  current  approach 
matches  those  from  Duntley’s  well  to  at  least  x  =  15. 

4.  Summary 

A  semi-analytical  point  spread  function  that  is  related  to  Wells/Gordon’s  formulation  is  presented 
and  show  close  fit  to  two  independent  measurement  results  under  different  conditions.  Monte  Carlo 
simulated  PSFs  are  calculated,  using  MVSM-measured  phase  functions.  Both  results  tested  well 
against  previous  empirical  results  from  Duntley  at  various  optical  lengths  and  scattering  contributions. 
The  fact  that  different  scattering  agents  can  result  to  an  order  of  magnitude  difference  in  PSFs  between 
measured  MC-MVSM  and  Duntley’s  at  small  angles,  or  in  other  words,  affecting  the  slope  of  the  PSF 
significantly,  suggest  additional  studies  are  needed,  particularly  in  natural  oceanic  waters. 
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